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INVESTIGATION  OF  THE  SPECTRAL  DECOMPOSITION  OF 
QUADRATURE  ERROR  SIGNALS 


1.  Introduction 

Two  aspects  of  radar  anti  communications  problems  are  the  de'ecti on  and  estimation 
of  received  signals.  Detection  is  'he  process  whereby  the  presence  of  a  received  signal  .s 
confirmed.  Estimation  seeks  to  quantify  the  detected  signal.  The  process  of  detection  and 
estimation  can  be  aided  by  use  of  a  signai  processing  technique  known  as  spectral  decom¬ 
position.  This  process  decomposes  an  observation  by  a  basis  ->et  which  spans  the  signal 
spaced  One  examDle  of  this  is  radar  doppler  processing,  whereby  the  presence  and  velocity 
of  a  target  can  be  determined.  When  the  particular  c'ass  of  signals  being  investigated  is 
periodic  the  choice  of  a  basis  set  for  the  decomposition  of  these  signals  may  lead  to  the 
selection  of  simple  periodic  functions,  the  sines  and  cosines. 

The  most  widely  used  technique  for  spectral  decomposition  based  upon  this  nasis  set 
is  the  Fourier  transform.  With  the  discrete  Fourier  iransform  (DFT)  a  set  of  N  uniformly 
spaced  samples  of  the  observed  signal  are  processed.  The  result  is  a  set  of  \  uniformly 
spaced  harmonic  estimates  of  the  signal's  periodic  spectrum.  In  practical  applications  of 
the  DFT  the  signals  of  interest  are  subjected  to  three  operations,  sampling,  the  DFT.  and 
windowing,  not  necessarily  in  that  order. 

The  discrete  Fourier  transform  has  been  examined  in  some  detail  and  some  pitfalls 
encountered  when  using  the  algorithm  have  been  investigated.  In  particular,  work  on  the 
problems  of  aliasing,  scalloping  loss,  and  spectral  leakage  can  be  found  in  the  literature.2 
This  paper  examines  the  special  case  of  spectral  decomposition  of  complex  signals  contain¬ 
ing  quadrature  error  perturbations.  Earlier  work  in  this  area  revealed  that  a  DFT  of  a 
complex  signai  vuh  quadrature  errors  would  generate  a  false  target  response  located  at  the 
mirror  image  position  relative  to  that  of  the  true  target.3  These  results  are  verified  and 
expanded  upon  in  this  paper.  Also,  the  following  algorithms  are  examined  for  possioie  use 
as  band-partitioning  filters  of  signals  with  quadrature  errors;  the  discrete  Hartley 
transform,  the  Hadamard  transform,  and  the  Real  DFT.  The  Real  DFT  is  the  ;erm  used  to 
describe  the  discrete  Fourier  transform  of  reai  data. 

2.  Problem  Definition 

Many  signal  processing  algorithms  reouire  digitized  nput  data  to  be  represented  bv 
complex  numoers.  These  complex  signals  can  oe  derived  as  part  or  'he  receiver  detection 
process  by  separating  real  signals  into  in-phase  mu  quadrature  components.  In  u  coherent 
ieiector  an  intermediate  frequency  signai  A  cos  '  to,,  -  nd!  .s  divided  equally  between 
two  oaseoand  mixers.  These  two  signals  ire  'hen  mxeu  w>th  •  multiplied  tv  •  he  'wo 
coherent  ocai  oscillator  tignais  given  ov  2  ..os >uU  ,mj  .1  mihu,.;  vspetuveiv.  The  -'ig- 
nais  exiting  from  tne  two  mixers  .re  'hen  low  pass  altered  \>  -eta  in  wuc  -tie  difference  :  re- 
ruenev  'fa,  component  oi  tne  mixing  process  Tie  -esuiiuK’  ignais  -ire  given  by 

/ (r  =  A  .os 
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These  signals  can  be  represented  as  a  complex  signal  he 
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Ideally,  or  detection,  the  t  wo  channels  a.  ili  idler  n  pnase  :  r<»m  one  anouier  ov  <u 
degrees  and  nave  the  same  amplitude  coefficient  in  reality,  the  pnase  and  amplitude  -an 
vary  greatly  depending  upon  the  quality  of  the  components  used  to  build  the  receiver  and 
on  environmental  conditions.  These  errors  are  attributable  to  RF  and  IF  filter  response 
mismatches.  RF  and  IF  timedeiavs.  quadrature  time  delays',  mixer  phase  and  amplitude 
mismatenes.  and  lowpass  tiller  mismatches,  to  name  a  few .t  For  this  investigation  ihe 
errors  will  be  modeled  as  phase  ind  amplitude  errors  .n  '.he  0  cr  innei  relative  to  the  1 
channel.  Phase  quadrature  error  is  the  amount  of  deviation  from  'he  .deal  pnase  separa¬ 
tion  of  90  degrees  between  the  I  and  Q  cnannels.  The  amplitude  quadrature  error  ;s  tne 
difference  in  gain  between  the  two  channels.  The  quadrature  signal  equations  can  now  be 
modified  to  reflect  these  errors  and  are  given  by 


/(. r  )  =  A  cos  Loxi 


Q(t  )  =  ( 1  -*-e  )A  sin(  to j£  —  6) 


wheie 

£  is  the  fractional  amplitude  error, 
O  is  the  phase  error. 


The  net  effect  of  these  phase  and  amplitude  quadrature  errors  on  a  signal's  spectrum 
is  to  create  a  false  target  response  at  the  frequency  which  is  the  negative  of  io\.  The  si/e  of 
this  error  response  has  been  quantified  by  Sinsky  and  Wang  and  a  technique  for  correcting 
these  quadrature  errors  has  been  proposed  by  Churchji],  Ogar,  and  Thompson5  and  ai.so  by 
Ifevdemann.0  The  power  of  the  false  image  relative  to  the  peak  of  the  ideal  -expense  is 
given  approximately  bv 
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image  power  (amplitude  error; 
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image  power  (phase  error) 


where  e  is  the  fractional  amplitude  ->rror  and  cf>  is  the  phase  error.  Churchill.  Ogar.  and 
Thompson  give  an  example  of  the  magnitude  of  the  effects  of  these  two  errors.  For  a  gain 
error  of  0.1  dB  and  a  phase  error  of  1  degree  the  image  power  wouid  appear  40  dB  below 
the  response  or  the  true  target.  In  practice  'hese  two  errors  can  approach  2  degrees  m 
phase  and  0.3  dB  in  amplitude.  Hrrors  of  this  order  can  cause  a  measuraole  degradation  in 
the  performance  of  many  signal  processing  algorithms. 

Once  a  detected  signal  oas  teen  ;pi it  into  1  and  Q  baseband  components  t  s  often 
iesiraoie  to  perform  some  vpe  >f  nana- partitioning.  The  discrete  Fourier  ransiorm  .s  a 
classic  'ecnmuue  for  resolving  l  sampled  waveform  ,nto  its  mectrai  components.  A 
oand -partitioning  filter  wmen  operates  on  real  .aiued  waxeiorms  is  known  as  'he  discrete 
Hartley  transform  12HT1.  knottier  ortnogonanzaiion  algorithm  'o  "e  investigated  is  a 
possible  'ransform  into  the  reuuency  domain  is  'he  Hauamard  transform.  The  liaiiatnaru 
transsorm  s  >i  interest  because  t  -an  -e  tnpiemenied  .-erv  <isnv.  :a c:i  >i  'tie  mentioned 
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The  inverse  DFT  is  given  by 

v_i  j2rryr 

/  (r)  =  £F {v)e  v  .  t-  0.1.2.  .V-l 

v-  0 

The  DHT  tor  real  valued  sequences  is  given  by 

.v  -1 

.FfiO  =  .V-1  £  /  (r )  caj  (277,i>t).  =  0.1.2.  .V— 1 
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The  inverse  DHT  is  given  as 
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The  Hadamard  transform  is  given  by 

F  ~  [H  }s  f 
f  =  N~1[H]nF 

where  N  -2k  k  =1,2 .  For  example 

t  H]2=  +  + 


[H]  4  = 


and 
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1+  ~  +  “ 
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It  is  interesting  to  note  that  both  the  DHT  and  the  Hadamard  transform  have  the  same 
weighting  terms  for  the  forward  and  inverse  directions.  In  practical  terms  this  couid  pro- 
videj.  significant  savings  m  hardware  since  the  same  weights  and  even  the  same  algorithm 
can  be  used  for  either  the  forward  or  inverse  transform. 


v  /  v  v.v 


3.  Effect  of  Quadrature  Errors  on  Band  Partitioning 

The  problem  s  lo  find  in  algorithm  which  ■'•ill  take  a  ..-mpiex  nru'  va. -d  mm  iiw 
iecemro.se  .t  info  is  spectrum  inu  el  be  tolerant  m  juaiiruture  *rr  >rs.  i-  air  ppu  umiv. 
nave  been  given  whicn  can  be  used  as  oand -partitioning  tillers.  It  remains  to  seen, 
'vhich  one  of  the  tour  can  be  used  in  a  system  that  has  F  Q  shanne!  perturbations  ^aused 
by  quadrature  phase  and  amplitude  errors  and  still  provide  a  desired  level  oi  performance. 
It  has  already  been  noted  that  the  complex  DFT  will  generate  a  false  response  which 
increases  as  the  square  of  the  magnitude  of  the  quadrature  errors.  This  phenomenon  will 
be  investigated  analytically  for  any  insight  it  might  provide  into  the  behavior  of  the  other 
transforms.  The  filters  will  also  be  examined  by  building  computer  models  >f  each 
transform  and  then  testing  the  band-partitioning  performance  of  each  algorithm  in  the 
presence  of  phase  and  amplitude  quadrature  errors. 

A  DFT  is  aiwavs  used  with  some  type  of  window  weighting  of  the  input  data.  For 
the  analysis  to  follow  a  rectangular  window  (uniform  weighting..1  is  used.  By  using  a  rec¬ 
tangular  window  the  frequency  domain  response  of  the  filler  will  be  realized  by  the 
coherent  addition  of  Dirichlet  kernels.’  The  first  input  waveform  to  bp  examined  is  a 
discrete  complex  sinusoid  given  by 


/  (n  )  =  e a 0  where  9  0  =  n  =  0. 1 .2.  .V  —  1 

The  discrete  Fourier  transform  of  this  complex  waveform  is 

F(9)  =  ,V-l'£Vn*0e--"3  = 

a  =0  n  =<) 


J  i‘v_1 

sin 

0) 

sin 

i-O  -e  0) 

2 

The  result  is  a  single  Dirichlet  kernel  centered  about  the  frequency  9  ,}  (see  Fig.  3).  This  is 
the  response  for  a  complex  sinusoidal  input  signal  if  no  quadrature  errors  are  present. 

A  discrete  complex  cosine  signal  with  phase  and  amplitude  quadrature  errors  can  be 
represented  as  follows 


I  (t  )  =  A  cos(n  9 
Q{t  )  —  (1  -Fe  ),4  sin(n  9  ,,  —  <p) 

!  f  t  )  -r  Q  1 1  >  =  A  cosi n  9  ,,)  +  j  (  1  t-  g  M  sin(  n  9  ,,  —  <i>  i 
where 


<S  s  the  fractional  amplitude  error. 


o  is  :he  phase  error  .*r  -am pie  Jeiav 

'Vith  each  error  "nken  •ndivuuaUv  and  lettim?  . ?  equal  unitv  ne  error  -*tiuation  '^c.-rnes 

'  n  '  =  cosvn  9.J  +■  ;  sim  n  tf  „  — -s.  quadrature  phase  error 
;  j  in  )  =  cost  n  0  .>)  *  i  ( 1  +€  )sim  n  9  quadrature  amplitude  error 


The  DFT  of  the  phase  quadrature  error  complex  signal  is  given  as  'details  m  the  Appendix! 
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The  discrete  Fourier  transform  has  revealed  an  input  signai  composed  of  four  fre¬ 
quency  components,  two  each  at  ±j9q.  The  two  kernels  attributable  to  the  imaginary  com¬ 
ponent  of  the  input  signal  have  a  scaling  factor  which  subjects  the  kernels  to  a  uniform 
twist.  This  result  was  obtained  by  using  the  shift  theorem  which  states  that  a  positive 
shift  of  a  function  in  time  will  result  in  a  uniform  twist  of  the  spectral  components  of  the 
signal  in  frequency.  Note  that  if  there  is  no  shifting  of  the  input  signal.  <f>=0.  then  the 
spectrum  reduces  to  that  derived  for  the  no  error  case.  In  fact,  the  response  for  the  no 
error  ca^e  is  composed  of  these  four  Kernels  which  coherently  combine  into  a  single 
response  at  9  0  wnen  o= 0. 

The  effective  spectrum  of  the  phase  uuadrature  error  signai  is  the  coherent  addition 
of  the  four  Diricniet  kernels.  It  can  be  seen  that  as  <i>  increases  the  ITeouency  response  of 
the  band-partitioning  filter  will  differ  markedly  :rom  the  ideal  response.  The  primarv 
response  centered  about  9.,  is  generated  by  the  addition  oi  two  kernels,  vhereas  the  mage 
response  centered  about  —9,,  is  generated  by  the  subtraction  ot  two  kernels.  One  of  each 
oi  .he  kernels  at  viil  ->e  twisted  ov  the  effect  af  the  ume  ,mft.  As  ’he  twisting 

.ncreases  the  kernels  at  —9,,  wtil  no  longer  ne  matched  and  he  ueep  null  at  this  :  reouencv 

will  oeg'.n  to  all  m  creating  an  image  response  jf  the  oruer  At  the  freauencv  -9  .  the 
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Similar  steps  are  taken  to  determine  the  DFT  ot  this  signal  as  were  done  lor  the  phase 
error  signal  ideiaiis  in  the  Appendix) 
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Again  the  four  Dirichlet  kernels  are  seen,  however  this  response  differs  from  that  of 
the  phase  error  signal  in  two  wavs.  First,  the  twisting  factor  due  to  the  lime  shift  is  no 
longer  present  and  second,  two  of  the  kernels  have  scaling  factors  other  than  unity.  As 
before,  if  the  error  perturbation  is  not  present  then  the  response  of  the  complex  signal 
coherently  combines  into  a  single  kernel  response.  The  net  effect  of  the  amplitude  error  is 
the  same  as  for  the  phase  error.  A  false  target  will  appear  at  the  mage  frequency  of  the 


true  target  with  a  magnitude  given  approximately  by 


isee  Fig.  15 ). 


ft  has  been  shown  that  a  complex  OFT  will  resoive  a  quadrature  error  signal  into  its 
SDectrai  components.  However,  in  addition  to  ’he  desired  main  response  there  will  be  an 
image  frequency  response  proportional  to  .ne  amount  it  quadrature  error  present  m  the 
channels.  The  image  frequency  response  s  caused  by  one  of  two  things;  a  twisting  of  the 
spectral  components  due  to  the  time  shift  .a used  oy  the  phase  error  or  an  amplitude  scal¬ 
ing  of  the  kernels  caused  by  the  amplitude  error,  in  a  .complex  DFT  -.he  conerent  addition 
of  the  kernels  takes  piace  within  the  transform.  To  avoid  the  comoining  it  the  quadrature 
error  terms  vumn  i  single  transform  a  r<* ai  DFT  is  used  on  eaon  ot  the  1  ana  Q  enanneis 
of  the  receiver.  The  iisauvantage  it  this  iDproack  ...  'hut  t.»e  .  - 1 ; n r» i  vnen  led 

into  two  reai  DFT’ s  will  transform  mo  our  -nit  mu  streams;  ;n  1  .hannei.  -eai  .nd  raa- 
gmarv  component  and  i  0  cnannei.  reai  inu  maginarv  comooneni  see  rig  a  The  advan- 
•age  cf  this  torroacn  s  nat  the  juadrature  jrror  ...mniex  lgnai  viii  be  resolved  nto 
jncomninea  -pectrai  terms.  That  as.  he  mauraiure  -rrors  viil  oe  soi.uen  » ithm  *acn  oi 
tne  'wo  f  and  D  enanneis. 


nr  examrte.  ;!  here  are  quadrature  Jrr  ms  r.  m  emit  :gna!  men  me 
v  lit  t>.;\  icr-'ss  me  1  .nu  <)  cunneis  l  me  .i rit : >  ;  in  S  .-■  'in>  ' )|  'r 

and  4  "  bands.  The  response  of  the  sec. mu  >m  >t  an  s  point  Dl  "  ;s  ;er 
nrut  samples  as  mi’,  ".vs 

k  -  Rt  =  i  ;  4-  <J  2  —  ‘  i  ~  C  j  -  l i.  ~  ‘  -  ~  C  i 
f2im  =  Qi-I2-Q,+  I^Qi-Jn~Q^U 

As  can  be  seen,  the  tutput  of  the  band  will  be  made  up  of  combinations  if  !  and  Q  sam¬ 
pled  data  points.  If  a  quadrature  amplitude  error  a  s  now  .ntroduced  mto  ’he  <J  channel 
the  response  of  the  imaginary  component  of  the  second  bin  will  be  : 

F  2|m  =  aQ  !  -  1 2  -  a<2  j  +  /  4  f  aQ  <  -  /„  -  aQ  -  t  l  s 

The  error  has  been  dispersed  unevenly  among  the  terms  and  therefore  cannot  be  easily 
removed  from  the  response.  A  possible  method  for  combating  this  problem  is  to  use  a  real 
DFT  on  each  of  the  I  and  Q  channels  of  the  receiver.  The  real  DFT  of  the  I  and  Q  channel 
data  is  given  as  follows  : 

F  2/r«  =  /  ;  —  /  3  +  /  5  —  /  -  :  r  71 1,.,  =  -  I  2  +  I  4~  ln  • h/s 

F2Q*'  =  Q2-Q<AQh-Q%;  F2Qlm  =  Q  ,  -  Q>  +  Q  5  -  Q  - 

The  presence  of  the  quadrature  amplitude  error  in  the  Q  channel  will  cause  me  real  and 
imaginary  components  to  be  scaled  as  shown  . 

F2Qi'  =  czhz-Q^Qn-Q*] 

F2Q{m  =  aj2,-Q,  +  Qs-C2r| 

By  transforming  each  l  and  Q  channel  individually  as  real  data  it  is  possible  to  isolate  the 
quadrature  error.  After  this  step,  an  appropriate  signal  processing  scheme  can  be  used  to 
compensate  for  these  errors.  A  more  general  discussion  of  the  DFT  of  real  data  follows. 

The  real  DFT  ot  the  discrete  complex  waveform  is  accomplished  n\  transforming 
each  channel  individually.  The  real  part  of  the  discrete  oompiex  input  waceform  is 

!  1  n  ,  —  cosin  Q  ,j.  n  -  0.1.2.  .V  —  1 

The  real  DFT  it  this  signal  :s 
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The  •  -  _trum  consists  of  two  kernels  centered  about  ±0  These  ire  the  cme  kernels  as 
1  r  obtained  with  the  complex  DPT.  The  difference  is  that  the  kernels  are  presented  as 
the  real  and  imaginary  parts  of  the  transform  of  the  1  channel  of  the  receiver. 

The  portion  of  the  phase  quadrature  error  waveform  in  the  Q  channel  of  the  receiver 
was  given  by 

f  ,  in  )  =  i  sinin  J.itiji 


The  real  J P  1  if  this  signal  is 
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Again  the  two  kernels  are  the  same  as  those  obtained  bv  using  the  complex  DPT.  8v  using 
the  real  DFT  and  band-partitioning  each  i  and  Q  channel  separately  the  resulting  kernels 
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4.  Computer  Results 


The  results  obtained  analytically  for  the  spectral  mp>s.;;,  n  i  i  rjadrar.re  armr 
siqr.ai  .an  :e  eru.eu  graphical!'.'  w  using  :ie  rmru'er  i"::e  '.mo  ir'  :i  icniri'  ’  lit-rs  .r-> 
examined  in  'he  following  order  the  complex  DFT.  he  reai  DFT.  .tie  DHT.  mu  he 
Hauamard  transform.  The  DHT  ana  iladamard  transform  are  being  .nvestigateu  because 
they  are  both  real  one  port  operations  which  can  be  built  to  run  faster  than  the  FFT.V  J  If 
they  have  satisfactory  band-partitioning  characteristics  they  can  be  used  instead  of  the 
real  DFT  for  decomposing  a  complex  signal  into  spectral  components. 


The  DFT  of  a  complex  valued  cosinusoidal  signal  should  produce  a  single  Dinchlet 
kernel  response  at  '.he  frequency  of  the  input  waveform,  in  Figure  i  3  <  the  response  oi  'he 
second  nm  \k  —2;  of  an  3  point  DFT  is  snown.  This  piot  is  obtained  bv  miecting  a  .,im- 

i  ’ 

piex  cosine  waveform  i  swept  in  frequency  from  —  2__  to  i  nto  the  S  point  complex 

/  . 

DFT.  For  the  values  selected,  bin  2  sr.ouid  have  its  peak  response  at  — —  and  mat  s  what 


:s  seen  :n  the  figure.  The  shape  of  the  kernel  is  determined  by  the  weight  applied  to  '.he 
sampled  data.  For  this  Investigation  a  rectangular  weight  is  used  and  this  is  confirmed  by 
the  -13  dB  sideiooes.  The  overlay  of  all  of  the  -S  bands  is  shown  :n  Figure  1  4  > .  This  plot 
shows  that  a  complex  DFT  will  decompose  an  input  waveform  into  S  filter  banks  each 


t , 

centered  about  the  frequence  k~~  where  /  is  the  sample  frequency  and  .V  is  the  number 

.V 

V  .V 

of  points  in  the  DFT.  The  index  k  runs  from  —  —  to  —  This  is  the  band-partitioning 


niter  response  for  complex  valued  waveforms. 


A  phase  quadrature  error  of  10  degrees  is  now  introduced  into  the  complex  input  sig¬ 
nal.  The  frequency  response  of  all  eight  bins  is  shown  in  Figures  (5-12).  The  frequency- 
response  of  the  second  bin  of  the  complex  DFT  is  shown  in  Figure  (7).  The  false  target 
response  is  cleariv  visible  at  the  image  frequency.  A  comparison  with  Figure  31  shows 
that  the  response  at  the  image  frequency  should  be  a  deep  null  and  not  a  -20  dB  response. 
The  reason  for  the  change  in  the  response  of  the  complex  DFT  is  that  the  coherent  addition 
of  the  four  Dinchlet  kernels  no  longer  results  in  a  single  kernel.  The  quadrature  error  of 
10  degrees  has  caused  the  spectral  components  of  the  kernel  to  be  twisted  uniformly  with 
frequency.  This  twisting  of  the  frequency  response  is  uue  to  the  time  delay  caused  bv  the 
phase  error.  Recall  that  at  +<# ,)  the  two  kernels  add  together  and  at  the  kernels  sub¬ 
tract.  The  result  of  this  twisting  of  the  kernels  is  most  noticeable  at  the  image  frequency. 
'Vhiie  there  may  not  be  a  noticeable  change  n  me  shape  of  the  response  at  the  primary 
frequency  :i  is  occurring 

The  amplitude  quadrature  error  will  produce  a  similar  change  in  the  complex  DFT’s 
frequency  response.  The  frequency  response  of  ail  the  bins  is  show  n  n  Figures  <  1  3 - 20 ' . 
For  these  riots  'he  amplitude  of  the  Q  channel  s  twice  that  of  the  I  cnannei.  The  DFT  ;.f 
the  impiitude  quadrature  error  signal  resuits  in  four  Dinchlet  kernels  centered  about  ex)  , 
where  two  or  me  kernels  have  scaling  factors  proportional  to  'he  amplitude  error.  The 
response  >f  the  ;econd  bin  of  the  DFT  s  snown  n  Figure  :  !5T  Again  tie  dise  arget 
response  ,s  seen  to  be  present  centered  ibout  the  mage  Tequencv.  The  -esuit  q 
-cnerenilv  aiding  the  four  mismatched  kernels  is  seen  'o  ilfect  the  overall  response  >i  he 
.liter  oank.  The  mange  in  'he  snare  of  he  tnamione  ma  -ideiones  of  'he  target  s  suite 
. : sable. 
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The  results  obtained  vith  '.”.e  computer  agree  .vith  t; 
nataeinatica.  ‘rsuiti  re',  emeu  be  necnanics  .*1  v  r.ji  \ca 
computer  riots  nrov  ide  visual  evidence  01  ‘die  net  result  m 
Diricnie;  seme  is.  "he  freuuencv  “esrv'r.se  u  're  -j:u  DFT 
-eai  Jala.'  w-.ii  re  examinee  next. 
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The  real  DFT  is  applied  lo  the  1  channel  of  the  receiver  and  the  output  of  bin  2  is 
shown  in  Figure  2 1 ) .  This  display  dilTers  from  that  seen  for  the  complex  DFT.  The  real 
DFT  produces  two  kernels  centered  about  sd instead  of  the  four  seen  with  the  complex 
DFT.  The  net  response  of  the  complex  DFT  was  a  single  kernel  centered  about  9 ,,  because 
of  the  way  in  which  the  four  kernels  coherently  combined  within  tne  t-ansform.  in  the 
real  DFT  there  are  no  opposing  kernels  centered  about  —9  ,.  there  is  only  one  kernel  as  ..an 
be  seen  in  Figure  t21).  The  overall  frequency  response  of  the  filter  bank  is  the  result  ol 
the  coherent  addition  of  these  two  kernels.  The  overlay  of  ali  the  tiller  bank  responses  is 
shown  in  Figure  i22).  It  can  be  seen  that  the  overall  response  maintains  the  symmetry 
seen  to  be  present  in  the  complex  DFT. 

The  transform  of  the  Q  channel  of  the  receiver  will  produce  plots  similar  to  those 
already  seen  for  the  I  channel.  The  effect  of  the  phase  Quadrature  error  wiil  -;e  to  rotate 
the  real  and  imaginary  vectors  about  their  axis.*0  The  magnitude  response  will  appear  to 
be  the  same  at  each  phase  error  step.  The  amplitude  quadrature  error  will  scale  the 
response  of  the  litter  by  an  amount  proportional  to  the  error. 

The  disadvantage  of  using  the  real  DFT  for  spectral  decomposition  is  that  it  creates  a 
two-in  four-out  type  of  network  .see  Fig  1 J.  However,  't  has  been  shown  that  this 
approach  will  decompose  a  quadrature  error  complex  valued  signal  into  its  periodic  spec¬ 
trum.  If  it  is  preferable  to  maintain  the  two-in  two-out  configuration  then  the  discrete 
Hartley  transform  and  Hadamard  transform  will  satisfy  this  requirement.  It  remains  to 
be  seen  if  these  two  algorithms  are  able  to  band-partition  a  time  domain  signal  as  well  as 
the  real  DFT. 


The  DHT  has  been  suggested  as  an  alternative  to  the  real  DFT  for  several  reasons. 
The  primary  reason  is  the  savings  in  compute  time  afforded  by  the  single  set  of  weignts 
needed  for  both  directions  of  the  transform.  The  same  is  true  for  the  Hadamard 
transform.  The  response  of  the  second  bin  of  the  discrete  Hartley  transform  for  the  cosine 
input  waveform  is  shown  :n  Figure  f23).  For  the  real  DFT  of  the  same  input  waveform 


the  mainlooe  was  centered  about  — j- 


For  the  DHT,  the  mainiobe  is  not  centered  about 


— p  ana  is  not  symmetrical,  "he  same  'vpe  of  benavior  is  exhibited  by  the  other  banas  is 

well.  The  overlay  of  ail  these  niters  is  shown  m  Figure  24).  The  DHT  may  be  useful  is 
a  mathematical  alternative  to  the  .-eai  DFT  but  it  joes  not  nave  much  merit  as  i  nand- 
partitioning  niter. 


The  'ast  filter  to  oe  examined  is  the  Hadamard  transform,  "his  ligoruhm  s  .erv 
easily  implemented  due  to  "he  simplicity  >x  the  weighting  terms.  As  vith  'he  DHT  the 
same  weights  are  used  for  the  forwara  and  inverse  transforms.  Hoe  ! lad:  .naru  weights  are. 
either  -1  or  -!.  The  DHT  ind  the  Hadamard  transform  iave  imiiar  vhavior.  In  tact. 

*.  =  4  'hen  the  two  transforms  are  denlicai. 
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with  that  for  •  he  DHT  as  seen  .n  Injure 


I  he  'derail  response  u  the  Hadamard 


ransiorm 


jure1,  2 b  i. 


comparison  v  ,tr.  '  ;ie  :oui  response  I  ’  ne  JUT  w  iil 


highlight  the  similarities  between  these  two  algorithms.  1‘nfortunately  ’he  1  iaJamarJ 
ansiorm  also  ex  .omits  'he  poor  band -partitioning  performance  seen  a.  :u:  '  ie  JUT 


5.  Conclusions 


This  investigation  was  undertaken  to  examine  the  behavior  of  four  band-partitioning 
filters.  The  complex  DFT  and  real  DFT  were  analyzed  mathematically  in  iruer  to  obtain  a 
better  understanding  of  the  mechanics  of  each  algorithm.  In  the  complex  DFT  tne  spectrai 
components  of  the  signal  are  coherently  combined  internal  to  the  transform.  This  process 
can  be  avoided  by  using  using  two  real  DFT's  to  filter  the  same  input  signal  and  thus  seep 
the  I  and  Q  channel  data  separated.  The  Hartley  and  Hadamard  transforms  were  investi¬ 
gated  by  developing  computer  models  of  each  algorithm  and  then  testing  each  filter. 


The  complex  DFT  has  been  widely  used  as  a  band-partitioning  filter  in  many  digital 
signal  processing  systems.  As  others  have  shown,  there  are  problems  with  using  this  filter 
in  receiver  systems  with  quadrature  errors.  The  most  significant  effect  of  I/Q  channel  qua¬ 
drature  error  is  to  cause  a  false  image  response  to  appear  at  a  frequency  which  is  the  nega¬ 
tive  of  the  frequency  of  the  target  response.  This  behavior  of  the  complex  DFT  was 
verified  both  mathematically  and  graphically.  It  is  because  of  this  perturbation  of  the  sig¬ 
nal  spectrum,  caused  by  the  mixing  of  the  quadrature  errors  between  the  real  and  ima¬ 
ginary  components  of  the  spectrum,  that  the  complex  DFT  makes  a  poor  band-partitioning 
filter  under  these  conditions. 


The  Hartley  and  Hadamard  transforms  operate  on  real  data.  Each  algorithm  uses  a 
single  set  of  weights  to  transform  between  the  two  domains  of  time  and  frequency.  By 
using  only  a  single  set  of  weights  less  memory  is  required  for  storage  when  these  filters  are 
implemented  in  hardware.  The  basis  set  for  the  Hartley  transform  is  given  bv 
cas  9  =cos0  +sin0  .  For  the  Hadamard  transform  the  basis  set  is  quite  simple,  the  values 
±1.  Even  though  each  of  these  transforms  can  be  efficiently  implemented  and  do  provide 
the  desired  quadrature  separation  they  do  not  make  good  band-partitioning  filters.  The 
reason  for  this  is  that  the  input  signals  of  interest  are  composed  of  sines  and  cosines  and 
the  Hadamard  transform  basis  set  is  not  made  up  of  these  waveforms. 


When  the  data  from  the  I  and  Q  channels  of  a  receiver  are  transformed,  using  the 
DFT.  as  two  real  numbers  instead  of  a  single  complex  number  then  the  desired  band- 
partitioning  filter  is  obtained.  The  difference  is  that  by  using  quadrature  separation,  chan¬ 
nel  mismatches  can  be  isolated  to  a  particular  spectrai  component.  By  doing  this,  succeed¬ 
ing  digital  signal  processing  techniques  can  be  used  to  compensate  for  these  channel 
mismatch  errors  and  thereby  achieve  a  desired  level  of  system  performance. 
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Fig.  1  —  Bandpartitioning  Filter  Configurations 
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Appendix 


The  DFT  of  a  signal  with  a  phase  quadrature  error  is  derived  as  follows 
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By  using  the  shift  theorem  the  following  can  be  shown 
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he  DFT  m  a  signal  vith  an  amplitude  quadrature  error  *s  .lerned  as  h.ulows 
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